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ABSTRACT 



The dynamics of the steady-state Sweet-Parker-type reconnection is 
analyzed in relativistic regime when energy density in the inflowing region is 
dominated by magnetic field. The structure of reconnection layer (its thickness, 
inflow and outflow velocities) depends on the ratio of two large dimensionless 
parameters of the problem - magnetization parameter a ^> 1 (the ratio of the 
magnetic to particle energy-densities in the inflowing region) and the Lundquist 
number S. The inflow velocity may be relativistic (for S < a) or non-relativistic 
(for S > a), while the outflowing plasma is moving always relativisticly. For 
extremely magnetized plasmas with o > S 2 , the inflow four- velocity becomes of 
the order of the Alfven four-velocity. 



1. Introduction 

Magnetic reconnection is widely recognized as a very important phenomenon in many 
laboratory and astrophysical plasmas (Biskamp 2000, Priest & Forbes 2000). It has been 
studied very extensively over the last 40 years, and a very significant progress has been 
made in understanding this process. However, historically, reconnection was of interest 
mostly to space physicists studying the Solar corona and the Earth's magnetosphere and 
to researchers in magnetic confinement fusion. In all these environments, plasma flows 
are non-relativistic and the Alfven velocity is usually much less than the speed of light 
(equivalently, magnetic energy density is much smaller than the particle rest mass energy 
density). Therefore it is not surprising that most of the progress on the subject has been 
made in non-relativistic regime. 
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Over the last decade however, it has been recognized that magnetic reconnection 
processes are also of great importance in high energy astrophysics, where dynamic behavior 
is often dominated by super-strong magnetic fields, with energy density B 2 / (8n) larger than 
the rest energy of the matter pc 2 + e. The best studied (but yet not completely understood) 
case is magnetized winds from pulsars. Models of pulsar magnetosphere (Goldreich & Julian 
1969, Arons & Scharlemann 1979, Ruderman & Satherland 1975) predict that near the light 
cylinder most of the spin-down luminosity of a pulsar should be in a form of Poynting flux. 
Other possible examples of relativistic strongly magnetized media include jets emanating 
from magnetized accretion disks around Galactic black holes and neutron stars as well 
as Active Galactic Nuclei (e.g., Beskin 1997, Lovelace et al. 2002), magnetosphere of 
magnetars (Thompson & Duncan 1996, Thompson et al. 2002) and Gamma Ray Bursters 
(Lyutikov k Blandford 2002). 

Dissipation of such super-strong magnetic fields may play an important role both 
for the global dynamics of the system and as a way to produce high energy emission. 
Magnetic reconnection has been proposed as the mechanism for acceleration of pulsar 
winds (Coroniti 1994, Kirk & Lyubarsky 2001) and GRB outflows (Spruit et al. 2001), 
and as a dissipation mechanism in AGN jets (Romanova & Lovelace 1992), Soft Gamma 
Ray Repeaters (Thompson & Duncan 1996), GRBs (Lyutikov & Blandford 2002, Spruit 
et al. 2001). In case of pulsar winds there are strong arguments that effective dissipation 
of magnetic field is, in fact, needed to account for the global dynamics of the Crab nebula 
(Kennel & Coroniti 1984; see also Michel 1994, Coroniti 1990, Melatos & Melrose 1996, 
Lyubarsky & Kirk 2001). 

This provides the motivation for studying magnetic reconnection in strongly relativistic 
plasmas (to be defined below). Despite of the growing interest in relativistic magnetic 
reconnection, very little theoretical (let alone experimental!) work has been done on the 
subject so far. We are aware of only one analytical discussion of relativistic reconnection 
(Blackman & Field 1994) and two recent numerical works on particle dynamics in relativistic 
reconnection layers (Zenitani & Hoshino 2001, Larrabee et al. 2002). 

In any theoretical analysis of magnetic reconnection, one first makes a number of 
approximations, e.g. incompressibility, two- dimensionality, the absence (or presence) of the 
axial magnetic field, etc. Then one formulates the set of MHD equations in a dimensionless 
form where the relative importance of various physical processes is represented by certain 
dimensionless parameters. After that one then tries to build a qualitative description of 
the reconnecting system by a small number of (also dimensionless) characteristic ratios. 
For example, in the simplest Sweet-Parker model of reconnection (e.g., Priest & Forbes 
2000) one assumes incompressibility, uniform and constant resistivity r/, and so on, and 
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one finds that the principal dimensionless parameter governing the system's behavior is the 
Lundquist number S = VaL/t] ^> 1, where Va is the upstream Alfven velocity and L is the 
size of the system. The other dimensionless plasma parameter, (3 - the ratio of thermal 
pressure to magnetic field energy density, turns out not to be as important, at least in the 
first approximation. Correspondingly S*™ 1 plays a role of the small parameter on which 
all further asymptotic expansions and boundary layer analysis are based. One then seeks 
to find out how the dimensionless characteristics of the system, such as the reconnection 
layer's aspect ratio 5/L and the ratio of the incoming velocity to the Alfven velocity scale 
with S as S — > oo. 

The generalization of reconnection to the relativistic case requires an introduction of 
one more (in addition to S) principal dimensionless parameter that should describe how far 
in relativistic regime we are. This parameter is a, the ratio of the magnetic field energy 
density to the plasma's rest mass energy density in the inflow region. We are now in position 
to define what we mean by relativistic reconnection: we consider a case when magnetic 
energy density in the flowing plasma dominates over particle energy density, er 3> 1. During 
reconnection magnetic energy is dissipated and transformed into plasma thermal energy and 
later into bulk motion. Since a 1, energy per baryon becomes much larger than m p c 2 and 
so we get relativisticly hot plasma in the reconnection layer. Expansion of plasma along the 
reconnection layer will produce relativistic bulk motions downstream. In relativistic regime, 
er 3> 1, the Alfven velocity becomes relativistic, Va = Ja/(1 + a)c ~ c, so that under 



certain conditions (to be determined later) the upstream flow is expected to be relativistic 
as well. 

Thus, in case of relativistic reconnection there are two very large dimensionless 
parameters, S and a. As we shall see below, one gets two different regimes depending 
on the ratio of these parameters: in one regime the incoming (the upstream) flow is 
ultra-relativistic, while in the other it is non-relativistic. 

In this paper we present a relativistic generalization of the simplest model of magnetic 
reconnection — the Sweet-Parker model — to strongly relativistic plasmas. This is a 
simple two-dimensional resistive MHD model, presented in Figure 1. We assume that the 
reconnection layer has a rectangular shape with a width L and thickness 5 <C L. The 
width L of the reconnection layer is determined by the global system size and thus, for the 
purposes of studying magnetic reconnection, is a fixed prescribed quantity. Also prescribed 
are the magnetic field strength and the baryon density and pressure in the ideal-MHD 
inflow region above and below the reconnection layer. In contrast, the thickness 5 of the 
reconnection region, as well as some other parameters such as the plasma inflow and outflow 
velocity, are not prescribed and need to be calculated as a part of the analysis. 
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We basically follow the steps of the Sweet-Parker analysis while taking into account 
relativistic effects, such as relativistic contraction and inertia of magnetic field. As the 
plasma enters the reconnection layer it slows down coming to a halt at stagnation point. 
At the same time the magnetic energy is dissipated and converted into internal energy of 
pair-rich plasma. In the out-flowing region the plasma is accelerated by the pressure gradient 
in x direction reaching some terminal relativistic velocity 7 out . We assume that energy losses 
are not important and the total energy of a fluid element (or at least a large fraction) stays 
within this fluid element, providing the corresponding amount of pressure support. Though 
generally radiative losses may be important, we expect that the reconnecting plasma will 
be optically thick to Thomson scattering after it's temperature becomes weakly relativistic, 
T > 20 keV. Above this temperature an efficient pair production process will start (e.g., 
Goodman 1986) inside the reconnection current layer trapping the radiation. The role 
of pair production in increasing the optical depth is an interesting question in itself and 
deserves further study (Thompson 1994). It lies, however, outside the scope of our paper; 
here we simply assume that the plasma is optically thick inside the reconnection layer and 
so the released magnetic field energy cannot leave the system and is therefore available for 
accelerating plasma downstream. 



2. Relativistic reconnection formulation 

The basic equations include the relativistic Ohm's law and the relativistic dynamics, 
Maxwell's and mass conservation equations (Lichnerowicz 1967): 

n = o, (i) 

= °> (2) 

(pu% = (3) 

where 

T ij = {w + b 2 + e 2 )«V + (P + ~ J ^)g ij - - eV (4) 

is the stress-energy tensor, w is the plasma proper enthalpy, p is proper plasma density and 
p is pressure, b 2 = bib 1 and e 2 = eVj are the plasma proper magnetic and electric energy 
density times 4n, p is pressure, u l = (7, 7/3) are the plasma four-velocity, Lorentz-factor and 
three-velocity, g 1 ^ is the metric tensor, b, t = ^rjij k iu^F kl are the four-vector of magnetic field, 
Levy-Chevita tensor and electro-magnetic field tensor and Cj = u^F^ is the four-vector of 
the electric field. 

The choice of the stress-energy tensor deserves some discussion. The stress-energy 
tensor (4) is a full electro-magnetic field plus matter tensor; this is not the relativistic 
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magneto-hydrodynamic (RMHD) stress-energy tensor. The reason is that in the frame 
work of RMHD it is assumed that one of the electro-magnetic invariants is not equal to 
and electro-magnetic stress energy tensor can be diagonalized. Equivalently, this implies 
that there is a reference frame where electric field is equal to 0. In case of resistive RMHD 
such frame may not exist, since, generally, there are resistive electric fields in the plasma 
rest-frame either along the null magnetic line (this violates then the B 2 — E 2 > condition 
of the ideal MHD), or field aligned electric fields (this violates E • B = condition). In the 
full relativistic approach the resistive electric field contributes to plasma energy density, 
energy fluxes and stresses. Nevertheless, we can still define the plasma rest frame by 
requiring that in that frame the electric fields e are only of resistive nature. We then obtain 
the stress-energy tensor (4). 

In what follows we will be using both the rest frame quantities (b, e, p, w, p*, p* for 
renormalized magnetic and electric fields, pressure enthalpy and charge and mass densities 
as well as laboratory quantities B, E, p e , p (pressure and enthalpy are defined only in the 
rest frame). 



2.1. Relativistic Ohm's law 

Relativistic Ohm's law is 

f = (j u )u j + -P k u k (5) 
rj 

where 77 is the plasma resistivity (e.g. Lichnerowicz 1967). In 3-D notations (Greek indexes 
= 1,2,3) this gives 

,- s/( . = a^-J.(B.« 

p z 777 

J = (j°-(j-/3))7 2 /3+^(E+(/3xB)) (6) 



which can be written 



The form (7) of the relativistic Ohm's law shows that for (j • /3) = and (E ■ j3) = the 
relativistic effects change the conductivity 

(8) 
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This can be understood if one notices that in the plasma rest frame the electric field 
e = E/i, and, since for (j • (3) = the current in the rest frame is j 



2.2. Main equations 

For stationary flow the energy and momentum flux conservation can be rewritten 
V a ( 1 2 (w + b 2 )(3)=0 

V p (y 2 wP a pP - p5 al3 ) = F^jfs = p e E a + (j x B) \ a (10) 
The Maxwell's equations then become 

(V x E) = 

divE = 47rp e (11) 

The equation of continuity is 

div 7/ 3p* =0 (12) 

The above equations plus the equation of state form a system of 15 equations for 15 
variables f3, E, B, j, p e , p,p. 

Following the Sweet-Parker model we neglect possible field-aligned electric fields and 
currents. Then the only non-vanishing components of electric field and current are E z and 
j z . Since the velocity lies in the x — y plane from eq. (6) it follows that the charge density 
p e = j° = 0. This is an important simplification since generally for relativistic plasma the 
charge density can not be neglected. 



2.3. Magnetization parameter 

We assume that far in the incoming region plasma is cold and strongly magnetically 
dominated. In the incoming region, well outside the reconnection layer, the restive electric 
field is vanishing and can be neglected. We introduce a frame-invariant ratio of the rest 
energy density of magnetic field and particles 

b 2 

a = ^>l (13) 
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This definition of a is equivalent to the ratio of Poynting to particle fluxes in the incoming 
region (c.f., Kennel & Coroniti 1984). For example, for pulsar winds initially (near the light 
cylinder) a ~ 10 3 — 10 6 (e.g., Arons & Scharlemann 1979). 

We have defined the cx-parameter in a frame-invariant way in terms of the rest-frame 
quantities. Alternative definition involves the laboratory frame fields and densities 

B? 

<Tl = — = 7in<7 (14) 
P 

It is straightforward to express the results in terms of 07 instead of a. 



3. Flow along the velocity separatrix 

In non-relativistic reconnection an important model problem is the flow of plasma 
along the velocity separatrix x = (e.g., Priest & Forbes 2000). In this case a simple 
equation relating magnetic field strength and inflow velocity can be obtained using only 
Ohm's law (7). Solving this equation for a given velocity (or vise versa for a given magnetic 
field) profile would allow us to find an example of velocity and magnetic distributions. 

The relativistic Ohm's law (7) along the velocity separatrix can be written as 

j{y) = j z (x = 0, y) = ^ (E z - (3 y B x ) . (15) 

For a more compact notation we introduce B(y) = B x [x = 0,y) > and also 
f3{y) = —f3 y (x = 0,y) > (here we changed the sign so that j3 is positive for convenience). 
Then, the above equation for Ohm's law can be rewritten as 

j(y) = j z (x = 0, y) = ^ [E z + (3{y)B{y)\ . (16) 



In a steady case Maxwell's equation gives V x E = 0, and hence, using d/dz = 0, 

E z (x, y) = const = E . (17) 
In the ideal region above the layer, Ohm's law gives 

E = -(3 in B in , (18) 
where B in = B(y > 5), and An = P(y > 5). 
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Using the first of the Maxwell equations (11) in laboratory frame, we find 

= - 7 (£ + (3B) , (19) 

or, using our expression for E, 

\B = (3 in -pB, (20) 
7 

where B = B / B in . This equation should be supplemented by the boundary conditions 
B(y = 0) = 0, /3(y = 0) = and B(y -> ±oo) = ±1, /% -> ±oo) = ±A n . Equation 
(20) relates the two functions — inflow velocity and magnetic field along the separatrix. 
For example, for a given f3(y) the two boundary conditions for the first-order ODE (20) 
determine the solution B(y) and put a constraint on the parameters 77, j3- m , and 5. 

For example, eq. (20) can be resolved for B(y) for a given (3{y): 

B(v) 8- f V P(y')-f{y')dy' r y f V ' f}(y")-f(y")dy" 

B(y) ee ^ = ^e" J 5 / 5 7^)^' • (21) 

-Din ^ J 

The boundary condition _B(y = 0) = is automatically satisfied, while the condition 
B(y = 5) = 1 serves as an eigenvalue problem for 5(r),fli n ). 

As an example, consider a case where the inflow four-velocity, u = 7/5, is a linear 
function of distance along the y axis: 

u= U4 if|y|<* (22) 
[ u in if |y| > 5 

By definition, <5 is the scale on which the inflow velocity changes from initial u = u- m to zero. 

Introducing dimensionless parameter Y = 5/rj (Y = c5/rj in dimensional units) and 
rescaling coordinate y by 5, y = y/5, the magnetic field is then determined by 

£(££1) = «*r-*J*fi f^nJTT^Jy (23) 
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Parameter Y here is an implicit function of u m given by the condition B(l) = B h 
Numerical solution Y(u- m ) is plotted in Figure 2. The corresponding magnetic field profiles 
for different values of u- in are plotted in Figure 3. 

Several important conclusions can be drawn from this exercise: (i) both the four 
velocity and magnetic field have the same typical scale ~ 5; (ii) both in the non-relativistic 
(■u in 1) and strongly relativistic (-u in 3> 1) regimes the ratio Y = cS/rj asymptotically 
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becomes inversely proportional to the four- velocity of the incoming flow, Y ~ l/w in , in 
agreement with the Sweet-Parker theory (in the non-relativistic case). Asymptotic scalings 
of Y(u{ n ) for u in <C 1 and u in ^> 1 are analyzed in Appendix A. (iii) in the case of strongly 
relativistic inflow, there is a thin sub-layer near the neutral point, y — 0, with the thickness 
^nr ~ S/u- m , where the flow becomes non-relativistic. 1 Outside of this sub-layer, the function 
B(y) approaches a universal shape in the limit u in — > oo. Inside the sub-layer the magnetic 
field becomes linear, 

B ~ --Bm, for y < 5/u in (24) 
V 

Typical magnetic field in the non-relativistic sub-layer is 

B nr ~ —B- m ~ % < An (25) 

T/ttin «tn 



while the typical current density is 

-Bin B; r 



?7 7in5 



(26) 



Here lies a qualitative difference between the non-relativistic and relativistic reconnection 
layer. In the non-relativistic Sweet-Parker theory the thickness of the layer is defined by the 
magnitude of current density at the center y = 0. In the case of relativistic inflow (7i n ^> 1), 
the current flowing in the bulk of the flow j ~ B in /5 is much stronger than typical current 
on the midplane. 

The example in this section is an illustration only, invoked in order to demonstrate 
a possible relation between 5 and f3- m . We use here only Ohm's law to establish such a 
connection; in reality, the velocity profile needs to be determined self-consistently by solving 
the entire 2D problem, including the equation of motion. 



4. Relativistic Sweet-Parker model 

In this section we use Ohm's law and the conservation laws for energy and particle flux 
and do not solve the momentum equation. All the estimates below are made up to the 
order of magnitude, or, more precisely, how they scale with two parameters, S and a. 



1 Because of the very large proportionality coefficient between Y and u ln in the strongly-relativistic case, 
this sub-layer develops only for very large (of order 10 2 and greater) values of u; n , as can be seen from 
Figure 3. 
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we find 



Estimating the current in the reconnection layer 

B TinAn-Bin 

Jz ~ y (27) 

AnTin ~ J- C (28) 

(we have restored the velocity of light here to make the right-hand side explicitly 
dimensionless) . Introducing relativistic Lundquist number 

S = — > 1 (29) 
V 

eq. (28) can be rewritten 

AnTin ~ ~ (30) 

This is the first basic equation of the model. 

The energy flux and the particle conservation give 

7inAn(l + <T)p* m L = 7outAutP*ut 5 

TinPniAn^ = ToutPout Aut<5 (31) 

In writing down these equations we have assumed that magnetic energy is fully spent on 
acceleration of baryons in the downstream flow. This assumes that all the electron-positron 
pairs created in the reconnection layer have annihilated. 

Then it follows that 

7in (l + a) = 7 out (32) 

and 



L S'AnTii 
P*nAn „ , _J 



PoutAut L 



1 + a)- (33) 



Equation (32) relates the velocity in the outflowing region to the inflow velocity 
and magnetization parameter. For small o 1 eq. (32) reproduces the familiar non- 
relativistic result that the outflowing velocity is of the order of the inflowing Alfven speed: 
Aut ~ V2a = B/c^Airp = Va/c. In the relativistic regime, a ^> 1, the outflowing velocity 
is always relativistic so that Aut ~ 1- It is remarkable, that in relativistic reconnection the 
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Lorentz factor of the outflowing plasma is much larger than the Lorentz factor of Alfven 
waves in the inflowing region 7^: 7 out ~ <77 in = 7^4 7m, where 7^ ~ y/a. 

Eqns. (32) and (33) are general relations for four unknown quantities: 5, 7 m , 7 out 
and p* n /Pout- 111 or der to resolve this system we need to use the momentum equations 
and find the structure of the flow inside the reconnection region, but this is a difficult 
task. An alternative way to proceed is to assume incompressibility of the plasma in its 
rest-frame p* n = p* ut (c.f., Blackman & Field 1994). Incompressibility may be due to, for 
example, a significant longitudinal field component B z . It is expected that the assumption 
of incompressibility will be justified as long as the inflowing plasma velocity is smaller that 
the fast magneto-sonic velocity, which in our case is similar to the Alfven velocity. 

Below we show that for a given value of S the inflowing velocity increases with the 
magnetization parameter a reaching at a ~ S 2 the Alfven velocity. For such strong 
magnetization the compressibility of plasma will become important. For any a <C S 2 the 
incompressibility is expected to be a good approximation. 

The assumption of incompressibility gives in relativistic limit a 3> 1 

/3 in 7 in L = 7 out 5 (34) 

and 

r 

An = <T- (35) 

(here we dropped A>ut because we expect an ultra-relativistic outflow with /3 out « 1.) 

Eq. (35) and the approximate Ohm's law (30) form a system of two equations for the 
thickness of the reconnection layer 5 and the inflow velocity An- There are two generic 
regimes (limiting cases) of relativistic reconnection: (i) relativistic inflow 7 in 1 and (ii) 
non- relativistic inflow j3 in <^ 1. 



4.1. Non-relativistic inflow An <^ 1 

Using eqns. (35) and (30) we find 

An ~ ~ ~ Y (36) 

Thus, 

* 1 «■ 1 (37) 



^in 



L cr Sp h 
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and finally 




where 7^ = \[2o is the Lorentz factor of the Alfven wave velocity in the incoming region. 
Since by assumption $ n <C 1 the non-relativistic inflow velocity is realized for a S. 

Thus, for a given a the inflow velocity is inversely proportional to the square root of 
the Lundquist number, similar to the classical non-relativistic Sweet-Parker model. 



4.2. Relativistic sub-alfvenic inflow 1 C 7 in < v2ff 
When (3- m rs 1 from eq. (35) it follows 

S 1 

(39) 



5 1 

r^j — 

L a 
and 

r 2 

s 



a 

7out ~ 



7in ~ I (40) 

For consistency we need cr ^> S (since we have always assumed relativistic motion 7 ^> 1). 
The ratio of the inflowing plasma Lorentz factor to the Alfven wave Lorentz factor is 

7in yJ°P 



(41) 

7a o 

Since the inflowing plasma should be sub-alfvenic for the incompressibility assumption to 
hold, it is required that 

a < 2S 2 (42) 
Thus the incompressible relativistic inflow case is applicable for 

S < o < 2S 2 (43) 



4.3. Relativistic alfvenic inflow 7 in ~ V2cr 

If a > 2S 2 the required inflow velocity becomes of the order of the Alfven velocity 
and the assumption of incompressibility should break down. Since the inflow velocity 
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cannot exceed the Alfven velocity for causal reasons, we assume that 7 in ~ 7^ = \f2a. The 
parameters of the reconnection layer then become 

7out ~ V2a^ 2 
5 1 



Pin (° 1 



P*out V 2 S 



> 1 (44) 



5. Bohm diffusion 

In the preceding section we have derived how the flow variables (first of all the inflow 
velocity) scale with two parameters of the model, S and a. For astrophysical applications 
it is often useful to have a qualitative estimate of the maximum possible reconnection rate 
in a given system. The maximum reconnection rate corresponds to the maximum value 
for resistivity and thus the minimal Lundquist number. This resistivity may be estimated 
using Bohm's arguments that the maximum diffusion coefficient in magnetized plasmas 
cannot be much larger than tlv, where is the Larmor radius and v is the typical velocity 
of electrons (of the order of the speed of light in our case). Thus, 

c 2 

77 ~ — (45) 

UJb 

and we find 

S ~ — (46) 

ym /m 

AlTin ~ ^ (48) 

Note, that in this /3 in approaches the velocity of light the reconnection layer becomes 

microscopically thin. One may expect that the fluid picture will become inapplicable at this 
point. 

For non-relativistic inflow velocity the assumption of Bohm diffusion gives 

5^ 



'r L L 



a 



An ~ (49) 
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while for relativistic inflow velocity 

~ 

7in 

7in ~ o"^r (50) 

which requires ar L /L > 1. In the case of relativistic inflow velocity and under the 
assumption of Bohm diffusion the thickness of the reconnection layer becomes smaller than 
the external gyro-radius. 

Equations (48) and the two limiting cases (49) and (50) give a useful "order-of- 
magnitude" estimates of the potential efficiency of reconnection in relativistic plasma. 



6. Discussion 

We have considered the dynamics of the relativistic Sweet-Parker reconnection under 
the assumption that the inflow region's energy density is dominated by magnetic field. We 
have found three generic regimes depending on the ratio of the magnetization parameter 
a to the Lundquist number S: (i) non-relativistic inflow velocity, cr < 5 1 ; (ii) relativistic 
sub-alfvenic inflow velocity, SCaC 2S 2 ; (iii) relativistic alfvenic inflow velocity, a > 2S 2 . 
For the first two regimes plasma flow may be assumed incompressible, while for the alfvenic 
inflow velocity compressibility is important. 

An apparent drawback of our approach is that we did not solve in a self-consistent 
way both the momentum and energy equations. This is a common flaw of many models of 
reconnection based on the Sweet-Parker approach. One can say that the role of the energy 
equation is to determine how compressible the plasma is. Conventionally, the simplest 
Sweet-Parker model does not include the energy balance equation with all its subtleties 
arising from possible radiation and conducting cooling effects; instead, it just replaces the 
energy balance equation with the incompressibility condition. In the absence of strong 
cooling or when a strong axial magnetic field component is present, the compressibility 
effects are indeed not important, at least in a rough, order-of-magnitude analysis. One 
then combines the incompressibility condition with the momentum equation and Ohm's 
law to arrive at the Sweet-Parker reconnection scaling. In our analysis we also assume 
incompressibity but then we use the energy conservation instead of the momentum equation; 
in the absence of energy losses these two approaches are, of course, equivalent and lead to 
the same results. 

In spite of the incompressibility assumption our approach represents a step forward 
in understanding relativistic reconnection as compared to the view of Blackman & Field 
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(1994). Using the same incompressibility assumption they were able to determine only 
the ratio of inflow and outflow velocities, while we find both these quantities separately, 
expressed in terms of external magnetization and the Lundquist number. We have also 
found a possible structure of the relativistic reconnection layer. 

For astrophysical applications, we were able to provide order-of-magnitude estimates of 
reconnection rates in relativistic plasma (eq. 48), which should be used instead of the often 
ad hoc assumptions of reconnection rates in strongly magnetized plasmas of pulsar winds 
and AGN jets. An important result is that under certain conditions the inflow velocity may 
become relativistic insuring very efficient dissipation of magnetic energy. 

We would like to thank Eric Priest, Vladimir Pariev, Eric Blackman, Chris Thompson 
and Alissa Nedossekina for comments on the manuscript. This research was supported by 
the NSF grant NSF-PHY99-07949. DU would like to thank CITA for hospitality during his 
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A. Asymptotic scaling of Y(u in ) 

In this appendix we investigate the asymptotic behavior of the function Y(u- in ) in the 
case of the 4- velocity field u(y) described by equation (22), for two limiting regimes, u- m — > 
and Wi n — > oo. 

We start by rewriting equation (22) in terms of rescaled coordinate y = y/S: 

a -M = Yl (^-0B). (Al) 

Now let us integrate this equation from to 1 and use the boundary conditions 
B(0) = and B(l) = I. We get 

i i 

B(l) - B(0) = 1 = Yfiin J >y(y)dy -Y J u{y)B{y)dy . (A2) 

o o 



This preprint was prepared with the AAS IATeX macros v4.0. 
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For the velocity field given by equation (22), i.e. u(y) = u in y, we then have 

1 i 

y- 1 = An / i{y)dy - / yB(y)dy . (A3) 



Using the relationship 7 = yl + u 2 , the integral in the first term on the right-hand 
side can be computed exactly, with the result 



- f 1 "i. n /\a_ u ^ I 

J i{y)dy = — J VTT^du = - — - — - + - Hu in + + <) . (A4) 
111 



In the strongly relativistic limit u- in — > 00 this expression can be expanded as 

]l(y)dy= 1 Y + \^n hl + 0{l). (A5) 



As for the second term, we note that in the strongly relativistic limit there is a thin 
non- relativistic boundary sub-layer of thickness 5 nr ~ S/ui n ; inside this sub-layer B(y) 
behaves linearly, with the slope that, as we shall see later, is inversely proportional to 
U[ n . However, what is important for us here, outside of this infinitesimally thin sub-layer 
the function B{y) approaches a certain universal shape, B TC \{y) in the limit u- m — > 00. 
This conclusion follows from our numerical solutions presented in Figure 3. Ignoring the 
non-relativistic sublayer's contribution to the integral in the second term on the right-hand 
side of equation (A3), we can then estimate this integral as 

1 1 

J yB(y)dy -> J yB rc i(y)dy, as u- m -> 00 , (A6) 


Combining this result with the result (A5) derived for the first term on the right-hand 
side of equation (A3), we finally can write the following expression: 

Y = - (A7) 

A re iu in + |lnw in + C" 

where we used the identity Jq 1 ydy =1/2 and defined 

1 

A-ei = j 'y[l- B rc i(y)]dy = const . (A8) 
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From our numerical solutions with very high (of order 10 3 ) values of u in , we find 
A Icl ~ 0.027 and C = -1. Thus, we get 

% 37 / 18.51nu in \ , A . 

Y(u- m ^oc)~— 1 ™ + ... . A9) 

M in \ M in J 

Notice that, because of the rather large (18.5) numerical coefficient multiplying the 
logarithmic factor, this logarithmic correction does not become negligible until one considers 
really very large values of u- m , of order 10 3 or more. This is the reason why the plot 
on Figure 2 (with u- m up to 200) does not show a very good agreement with the simple 
power law Y ~ l/u [n . Notice however, that the numerical results for values of u- in in the 
range between 200 and 2000 are in excellent agreement with the predicted asymptotic 
behavior (A7), as demonstrated separately in Figure 4. 

Now let us consider the non-relativistic regime, u- m < 1. In this limit, the 
expression (A4) for the integral entering the first term on the rhs of equation (A3) tends 
to 1 (this is of course a trivial result since in this limit 7(2/) ~ 1). Then, since f3 in ~ u in in 
this regime, this first term can be evaluated simply as u in . 

At the same time, just as in the strongly relativistic case, the function B{y) also 
approaches a certain universal profile, as can be seen in Figure 3 (we shall call this 
profile B nr (y)). Then, the integral entering the second term on the rhs of equation (A3) 
asymptotically approaches a constant value 

} 1 

J yB(y)dy -> - - A ni = const , (A10) 


where we defined A nr in a manner similar to A rc \. 

1 

A nT = J y[l-B m {y)]dy. (All) 


Thus, we see that in the non-relativistic limit the function Y(u in ) asymptotically 
approaches an inverse power law: 

1 1 1.71 
Y(u m ) * - T — r * — , (A12) 

where we used the value of A nT = 0.0855 computed using our numerical solution for the 
case -u in = 0.01. The resulting asymptotic behavior in the non-relativistic regime agrees 
very well with the numerical results, as can be seen in Figure 2. 
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Fig. 1. — Geometry of the model. The reconnection layer has a width L in x direction and 
thickness 5 in y direction (and infinite depth in z). Outside there is strongly magnetized 
plasma flowing into the reconnection layer with velocity along the y direction. Inside 




Fig. 2. — Dependence of Y = 5/rj on the four velocity of the incoming plasma for a linear 
dependence of u(y) = u- in y/8 inside the reconnection layer. The non-relativistic asymptotic 
isY = 1.7/iiin. 
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Fig. 4. — Asymptotic behavior of the function Y(u- m ) in the strongly-relativistic limit 
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